Measurement of cosmic-ray antiproton spectrum at solar minimum with a 
long-duration balloon flight in Antarctica 
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The energy spectrum of cosmic-ray ant.iprot.ons (p’s) has been measured in the range 0.17 to 3.5 
GeV, based on 7886 p’s collected by the BESS-Polar II instrument during a long-duration flight over 
Antarctica in the solar minimum period of December 2007 through January 2008. The p spectrum 
measured by BESS-Polar II shows good consistency with secondary p calculations. Cosmologically 
primary p’s have been searched for by comparing the observed and calculated p spectra. The BESS- 
Polar II result shows no evidence of primary p’s originating from the evaporation of PBH. 


The precise measurement of the spectrum of cosmic- 
ray antiprotons (p’s) is crucially important to investiga- 
tions of conditions in the early universe and of cosmic- 
ray propagation. Most cosmic-ray p’s are produced by 
energetic nuclear cosmic rays interacting with the inter- 
stellar gas. The energy spectrum of such “secondary” 
p’s is expected to show a characteristic peak around 2 
GeV. with sharp decreases of the flux below and above 
the peak due to the kinematics of p production and the 
interstellar proton spectrum. The secondary p’s offer a 
unique probe [1-3] of cosmic-ray propagation and of so- 
lar modulation. Other possible sources of cosmic-rav p’s 
have been sought, such as primordial black holes (PBH) 
evaporating by Hawking radiation [4], The “primary” 
spectral contributions of PBH sources, if they exist, arc 
expected to be observable at low energies [5] and to ex- 
hibit large solar modulation effects [6] due to the shape 
of the predicted spectrum. 

In the previous solar minimum period, the BESS ex- 
periment (BESS95+97) showed that the p spectrum has 
a distinct peak around 2 GeV, as expected, and it has 
become evident that p’s are predominantly of secondary 
origin [7] . However, the low-energy component of the p 
spectrum measured by BESS95+97 was slightly flatter 
than predicted by calculations with secondary p models. 
Though this result might suggest the existence of novel 
processes for production of cosmic-ray p’s, the large sta- 
tistical error of the BESS95+97 data did not allow a firm 
conclusion. The BESS-Polar project [8-10] was proposed 
to evaluate the possibility of excess low-energy p flux sug- 
gested in the BESS95+97 observation by measuring p’s 
with unprecedented precision using long-duration flights 
over Antarctica at solar minimum. The flight of BESS- 


Polar I was performed in December 2004 [11-14] and the 
flight of BESS-Polar II [12] was carried out near solar 
minimum during December 2007 and January 2008. In 
this Letter, we report a new measurement of cosmic-ray 
p’s in the energy range 0.17 GeV to 3.5 GeV by BESS- 
Polar II. Based on the measured spectrum, the subjects 
of secondary p’s and possible primary p’s are discussed. 

The BESS-Polar instrument [15-17] was designed and 
developed as a high-resolution magnetic-rigidity spec- 
trometer. A uniform field of 0.8 T is produced by a thin 
superconducting solenoid, and the field region is filled 
with drift-chamber tracking detectors. Tracking is per- 
formed by fitting up to 52 hit points with a characteris- 
tic resolution of ~ 140pm in the bending plane, result- 
ing in a magnetic-rigidity (= Pc./Ze) resolution of 0.4% 
at 1 GV and a maximum detectable rigidity (MDR) of 
270 GV. Upper and lower scintillator liodoscopes provide 
time-of- flight (TOP) and dE/dx measurements and the 
event trigger. For antiproton measurements, the accep- 
tance of BESS-Polar is 0.23 m 2 sr. The timing resolution 
of the TOF system is 120 ps, giving a 8~ 1 resolution of 
2.5%. The instrument also incorporates a threshold-type 
Cherenkov counter using a silica aerogel radiator with 
index n = 1.03 (ACC) that can reject e~ and p back- 
grounds by a factor of 6100 and distinguish p’s from such 
backgrounds up to 3.5 GeV. In addition, a thin scintilla- 
tor middle-TOF (MTOF) is installed on the lower surface 
of the solenoid bore to detect low-energy particles which 
cannot penetrate the magnet wall. The timing resolu- 
tion using the MTOF is 320 ps. In the present analysis, 
the MTOF was employed to verify the procedure used to 
eliminate contamination from events in which interacting 
protons mimic low-energy p’s. 
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TABLE I. p flux at the top of atmosphere with statistical (first) and systematic (second) errors. Np and Nbg are the number 
of observed antiprotons and estimated background events. 
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The BESS-Polar II payload was launched on December 
23, 2007, from Williams Field near the US McMurdo Sta- 
tion in Antarctica and circulated around the South Pole 
for 24.5 days of observation with the magnet energized. 
The float altitude was 34 km to 38 km (residual air of 5.8 
g/cm 2 on average), and the cutoff rigidity was below 0.5 
GV. BESS-Polar II accumulated 4.7 x 10 9 events with no 
inflight event selection as 13.6 terabytes of data. 

During the flight, most detectors and instrument sys- 
tems operated well and provided their expected perfor- 
mance. The central tracker exhibited high-voltage flue- 
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FIG. 1. The 0 1 versus rigidity plot and p’s selection band. 
The solid curves define the p mass bands. The lowest energy 
p’s are shown in the closeup figure. 


tuations. However, normal tracking resolution was ob- 
tained for more than 90% of the science observation time 
by development of algorithms that calibrate the tracker 
over short time intervals and depend on its high-voltage 
state. One upper TOF PMT and one lower TOF PMT 
had high-voltage control problems and were turned off. 
A conservative requirement of two good PMTs on each 
TOF scintillator resulted in a ~ 20% reduction in accep- 
tance. 

Analysis was performed in the same way as described 
in Ref. [7 . The same selection criteria for p’s and protons 
were applied because noninteracting p’s behave similarly 
to protons in the symmetrical configuration of BESS- 
Polar, except for deflection direction. 

Figure 1 shows /? -1 versus rigidity plots for surviving 
events. We see a clean narrow band of 7886 p’s at the ex- 
act mirror position of the protons. The ACC veto helped 
identify p events at higher energies by removing the back- 
ground of e~ and p~. The fractions of contamination 
were 0.0%, 1.0%, and 2.3% in the 0.2-1. 0 GeV, 1. 0-2.0 
GeV, and 2. 0-3. 5 GeV energy" bands. Other backgrounds 
such as albedo, mismeasured positive-rigidity particles, 
and re-entrant albedo were found to be negligible. 

After p candidates are identified, the differential flux 
of p’s at the top of atmosphere (4>toa) integrated in an 
energy" -width of d£ can be expressed as follows: 

■J'TOAdE = (A T TOI ~ Aatmos)/£air/(SD ' Tiive) (1) 
A'tOI = (Np — A/BG)/(edet ' £non-int) (2) 

where Ti; vc is the live time, Np and Affic are the num- 
bers of observed p candidates and expected background 
particles among the candidates. For the present anal- 
ysis Tiive totaled 1286460 seconds. The effective geo- 
metrical acceptance including non-interaction efficiency 
(5'n-e non _ int ) was calculated using GEANT3 to be 0.133 
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FIG. 2. Antiproton flux at the top of the atmosphere obtained 
with BESS-Polar II at solar minimum together with BESS 
flights around the previous solar minimum (BESS95+97) and 
PAMELA. 

± 0.011 m 2 sr at 0.2 GeV and 0.159 ± 0.008 m 2 sr at 2.0 
GeV. with errors estimated from differences relative to 
GEANT4 calculations. The detection efficiency of p's 
(sdct) w as calculated using a noninteracting proton sam- 
ple to be 81.4 ± 0.1 % at 0.2 GeV and 60.0 ± 0.2 % at 
2.0 GeV. Estimates of atmospheric secondaries included 
in A’toi (Aatmos) were calculated with the proton and 
helium fluxes measured by BESS-Polar II as input. The 
subtraction amounts to 17.6 ± 2.0 % at 0.2 GeV and 

27.6 ± 0.1 % at 2.0 GeV, where the errors correspond to 
the maximum difference among three recent calculations 
[18-20], In order to obtain the TOA flux, a correction 
for survival probability [19] in the residual atmosphere 
(e air ) was applied, and the probability was estimated as 

85.6 ± 2,0 % at 0.2 GeV and 89.8 ± 2.0 % at 2.0 GeV. 

Table I gives the resultant flux of p’s at the top of the 

atmosphere in the kinetic energy range 0.17 to 3.5 GeV 
with the statistical (first) and systematic (second) er- 
rors. The dominant systematics come from atmospheric 
subtraction and detection efficiency. A rapid change in 
efficiency due to the ACC veto increases the systematic 
uncertainty in the two highest bins. 

Shown in Fig. 2 is the BESS-Polar II p spec- 
trum together with measurements by BESS95+97 and 
PAMELA, and various secondary p calculations [18, 21- 
24]. Improved statistical precision of the p flux measure- 
ment results from 14 ahd 30 times more events below 1 
GeV than -were measured in BESS95+97 and PAMELA 
[25]. respectively. The p spectrum measured by BESS- 
Polar II generally shows good consistency with secondary 
p spectra calculated for solar minimum conditions. 

The evident differences among the calculations shown 


FIG. 3. Comparison of antiproton flux shapes with secondary 
p calculations. In order to focus on the difference of spectral 
shapes in the low-energy region, the calculations are normal- 
ized at the spectral peak ~ 2 GeV. The uncertainty due to 
the modulation parametrization with the range from 500 MV 
to 700 MV in the Mitsui model is shown by the shaded band. 

in Fig. 2 arise from several factors: (1) definition of 
the primary proton and helium spectra, (2) incomplete 
knowledge of nuclear physics in propagation, (3) parame- 
ters and models of propagation in the Galaxy, and (4) the 
modulation in the heliosphere. These factors affect both 
the normalization and shape of a calculated p spectrum 
in different ways and degrees. Variation in the absolute 
fluxes of interstellar protons and helium will, for instance, 
affect the absolute flux of p’s produced, but not the shape 
of the spectrum. 

In Fig. 3, we normalize all of the p spectrum calcu- 
lations near the peak energy at 2 GeV to focus on their 
shapes. The calculated spectra and data points arc also 
multiplied by Ep 1 to further emphasize differences at low 
energies. The data have not been normalized. Precise 
measurement of the p spectrum to low energy by BESS- 
Polar II allows secondary p calculations to be evaluated 
by comparing predicted and observed spectral shapes. 
Chi-square (x 2 ) calculated with BESS-Polar II data and 
normalized secondary p calculations in Fig. 3 are 0.57 
(curve 1), 0.56 (curve 2), 1.24 (curve 3), 1.59 (curve 4), 
0.63 (curve 5). The shaded band in Fig. 3, calculated 
using the Mitsui model [6] modulated with a standard 
symmetric approach (Force-Field), indicates the small 
variation that results from uncertainty in the modula- 
tion parameter. Shown is the range from 500 MV (lower 
border: y 2 =0.77) to 700 MV (upper border: x 2= 0.48). 
Even though Force-Field modulation over this range of 
parameters gives a visible change to the spectral shape, 
as does the Drift model, this is small in comparison to 
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FIG. 4. (Top) Possible primary antiproton fluxes from 
evaporating PBHs calculated for BESS-Polar II (A) and 
BESS95+97 (B) by fitting differences of the measured spectra 
from the Mitsui secondary p spectrum. (Bottom) Distribution 
of the evaporation rate of PBH (77) calculated for BESS-Polar 
II and BESS95+97. Negative values of 77 are non-physical. 


variations among the propagation models, because of the 
shape of the low-energv p spectrum. BESS-Polar II re- 
sults show better consistency with the models that do 
not include a component of low-energy p’s from sources 
such as tertiary p’s or a soft spectrum such as that due 
to the Diffusive Reacceleration model (curves 3 and 4 in 
Fig. 3). 

Among the various possible sources of primary p’s, 
the most pertinent candidate for BESS-Polar to study 
is the evaporation of PBHs that may have formed in 
the early Universe via initial density fluctuations, phase 
transitions, or the collapse of cosmic strings. PBHs are 
the only black holes that have small enough mass to 
emit particles with a significant “evaporation” rate at 
the current age of the Universe. The possible existence 
of primary p’s has been quantitatively evaluated by fo- 
cusing on the PBH evaporation rate parameter (77) cal- 
culated by fitting the differences of p spectra measured 
in BESS-Polar II and BESS95+97 from theoretical sec- 
ondary p calculations with normalization near the peak 
at 2.0 GeV. The normalization is performed to avoid 
bias from ambiguity in the absolute fluxes predicted by 
propagation models. The dependence on model choice is 


also considered by checking the primary p’s under sev- 
eral models. The BESS-Polar II spectrum and the con- 
servative models and modulations shown in Fig. 4 give 
1Z = 5.2^' j x 10~ 4 pe“ 3 yr _1 . This excludes by 8 sigma 
the slight possibility of primary p’s suggested by the 7 7 
distribution with a mean value of 4.4 x 10~ 3 pc'" 3 yr~ 1 
determined from BESS95+97 data. We have also deter- 
mined an upper limit of 77. ~ 1.2 x 10 " 3 pc~ 3 yr _1 with a 
90% confidence level. This upper limit is almost insensi- 
tive to the modulation used in the secondary model calcu- 
lation (500 MV: 77 = 1.1 x 10~ 3 pc _3 yr - \ 600 MV: 77 = 
1.2 x 10 _3 pc _3 yr _1 , 700 MV: 77 = 1.3 x lO^pc-V" 1 ). 
Within statistics, the BESS-Polar II result shows no evi- 
dence of primary p’s originating from the evaporation of 
PBH. 
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